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Abstract 

A  stability  analyzer  for  testing  NASA  Deep  Space  Network  installations  during  flight  radio 
science  experiments  is  described.  The  stability  analyzer  provides  realtime  measurements  of  signal 
properties  of  general  experimental  interest:  power,  phase,  and  amplitude  spectra;  Allan  deviation; 
and  time  series  of  amplitude,  phase  shift,  and  differential  phase  shift.  Input  ports  are  provided 
for  up  to  four  100  MHz  frequency  standards  and  eight  baseband  analog  (>I06  kHz  bandwidth) 
signals.  Test  results  indicate  the  following  upper  bounds  to  noise  floors  when  operating  on  100  MHz 
signals:  —145  dBcIHz  for  phase  noise  spectrum  further  than  200  Hz  from  corrier,  2.5  x  10 
(t  —I  second)  and  1.5  x  10“^^  (t  =1000  seconds)  for  Allan  deviation,  and  1  x  10  degrees  for 
1 -second  averages  of  phase  deviation.  Four  copies  of  the  stability  analyzer  have  been  produced, 
plus  one  transportable  unit  for  use  at  non-NASA  observatories. 


Introduction 

The  Deep  Space  Network  (DSN)  is  called  upon  to  attain  high  levels  of  frequency  stability  for 
scientific  purposes.  For  instance,  the  upcoming  Cassini  mission  to  Saturn  will  use  the  DSN 
to  attempt  detection  of  gravitational  radiation,  and  to  observe  properties  of  Saturn’s  rings, 
atmosphere,  and  satellites 91. 

These  and  related  investigations^!  measure  small  perturbations  on  a  radio  signal  passing 
between  the  earth  and  a  distant  spacecraft.  The  Cassini  applications  are  fairly  typical,  requiring 
frequency  stability  of  a  few  parts  in  10^'^  (Allan  deviation  for  sampling  time  r  =100  to  10,000  s) 
and  single-sided  phase  noise  around  -60  dBc/Hz  (1  to  10  kHz  offset  from  an  8.4  GHz  carrier). 

It  is  challenging  to  achieve  such  stabilities  in  the  operational  environment  faced  by  the  DSN. 
That  environment  includes  months-long  periods  of  duty;  spatially  distributed,  outdoor,  and 
moving  equipment;  and  competition  for  observing  time.  We  have  found  that  stability  failures 
can  remain  hidden  in  the  bulk  of  DSN  activities,  only  to  surface  when  the  scientific  experiment 
is  undertaken.  This  is  troublesome  because  most  mission  experiments  cannot  be  repeated. 
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Therefore  the  DSN  has,  in  the  past,  tested  its  systems  using  instrumentation  suitable  for  use 
by  specially  trained  personnel.  This  approach  was  expensive,  however,  and  the  time  to  analyze 
data  has  often  allowed  additional  diagnostic  evidence  to  disappear,  necessitating  repeated  tests. 

We  developed  a  stability  analyzer  to  enable  operations  personnel  to  rapidly  measure  stability  in 
various  ways,  in  order  to  lower  costs  and  reduce  response  time.  The  particular  measurements 
made  are;  power,  phase,  and  amplitude  spectra;  Allan  deviation;  and  time  series  of  amplitude, 
phase  shift,  and  differential  phase  shift.  Our  analyzer  provides  inputs  for  up  to  four  100 
MHz  frequency  standards  and  eight  baseband  analog  (>100  kHz  bandwidth)  signals,  with  the 
possibility  of  expanding  to  accept  digital  inputs  over  a  local  area  network.  Four  copies  of 
the  stability  analyzer  have  been  produced,  plus  one  transportable  unit  for  use  at  non-NASA 
observatories. 


Instrument  Overview 

TTie  DSN  stability  analyzer  has  two  major  components:  1)  the  RF  and  Analog  Assembly,  and 
2)  the  Cxintroller  Assembly,  as  depicted  in  Figure  1. 

The  RF  and  Analog  Assembly  provides  the  conditioning  and  conversion  of  the  input  analog 
signals  into  a  signal  the  controller  can  analyze.  The  equipment  is  installed  in  two  parts:  the 
100  MHz  Interface  Assembly  and  an  RF  Cabinet  Assembly. 

The  100  MHz  Interface  assembly  resides  as  close  as  possible  to  the  DSN  primary  frequency 
standards,  usually  hydrogen  masers  (H-masers).  Intentionally,  this  location  is  isolated  from 
routine  personnel  access,  as  well  as  from  as  many  environmental  influences  as  possible.  The 
assembly  receives  four  100  MHz  inputs,  which  are  compared  in  pairs.  The  comparison  (described 
further  below)  results  in  a  100  kHz  signal  that  is  sent  over  a  fiber-optic  interface  to  the  RF 
cabinet.  The  RF  assembly  resides  in  a  convenient  location  for  access  by  test  personnel.  It 
provides  reference  frequency  synthesis  and  distribution,  switching  among  the  possible  input 
sources,  signal  conditioning  in  the  form  of  amplification,  and  optional  downconversion  with 
detection  of  zero  crossings. 

The  Controller  Assembly  resides  next  to  the  RF  assembly,  and  provides  an  operator  interface 
for  selection  of  the  test  type  and  hardware  configuration,  and  for  presentation  of  results.  The 
Controller  also  controls  details  of  switches  and  instrumentation,  acquires  data  by  means  of 
analog-to-digital  (A-D)  converters  and  a  time  interval  counter,  and  analyzes  the  data  acquired. 
Originally,  the  RF  assembly  was  housed  in  one  rack  and  the  controller  equipment  was  housed 
in  a  second  rack.  These  cabinets  have  since  been  bolted  together  to  form  a  double  cabinet, 
and  components  of  each  have  been  swapped  to  improve  ergonomics  for  the  operator.  See 
Figure  2  for  a  photograph  of  the  double  cabinet. 
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Analog  Electronics  Design 

100  MHz  Interface  Assembly 

This  assembly  selects  the  pair  of  100  MHz  signals  to  be  analyzed,  and  converts  the  selected 
signals  into  a  form  that  can  be  transported  to  the  low  frequency  equipment.  Figure  3  shows  a 
block  diagram. 

Output  from  the  100  MHz  assembly  is  sent  on  fiber  optics  to  the  low  frequency  equipment 
to  prevent  ground  loop  currents  that  could  induce  spurious  signals  or  noise  into  signals  being 
mea.siired,  or  could  contaminate  the  frequency  standard’s  outputs.  The  100  MHz  Interface 
has  four  100  MHz  input  ports.  Two  input  ports  are  connected  to  H-Maser  outputs,  and 
one  other  port  is  normally  used  for  comparing  the  station’s  coherent  reference  generator  100 
MHz  output  against  the  H-masers.  The  100  MHz  signals  are  selected  for  measurement  using 
RF  relays  followed  by  high  reverse  isolation  amplifiers  cascaded  with  output  matrix  switches. 
The  combined  isolation  of  both  sets  of  switches  and  60  dB  reverse  isolation  of  the  amplifiers 
provides  more  than  150  dB  crosstalk  isolation  between  signals. 

Switch  control  commands  are  sent  over  fiber  optics  to  the  100  MHz  Interface  using  commercial 
modems  and  digital  I/O  boards  to  address  switch  decoders  that  operate  the  switches. 

The  selected  pair  of  inputs  are  frequency  multiplied  by  99  and  100  respectively  with  phase-locked 
cavity  multipliers.  The  multiplier  outputs  at  9.9  GHz  and  10.0  GHz  are  mixed  to  generate  100 
MHz.  The  result  is  frequency  translated  to  100  kHz  in  an  offset  frequency  generator,  and  sent 
on  a  fiber-optic  link  to  the  low  frequency  assembly. 

The  frequency  conversion  process  yields  a  single  100  kHz  carrier  with  a  phase  spectrum 
containing  the  relative  stability  of  the  100  MHz  inputs,  with  a  40  dB  margin  above  what  would 
be  obtained  from  direct  mixing  of  one  input  with  the  other  input,  offset  by  100  kHz.  Amplitude 
information  is  lost.  Frequency  translation  to  100  kHz  is  necessary  for  the  A-D  converter,  and 
allows  the  signal  to  be  transported  to  the  low  frequency  assembly  over  low-cost  multimode 
fiber  optics. 

Low  Frequency  Interface  Assembly 

This  assembly  contains  switches  that  select  among  baseband  receiver  signals  and  the  100  kHz 
.signal  from  the  100  MHz  assembly.  The  selected  signals  are  routed  to  measurement  ports  of 
the  computer  system.  Figure  4  shows  a  block  diagram. 

Baseband  signals  are  selected  by  matrix  switches  and  sent  to  programmable  attenuators  that 
set  levels  into  the  interface  amplifiers.  Another  matrix  switch  outputs  the  selected  signals  to 
the  desired  output  ports.  The  frequency  translated  lOO  MHz  maser-pair  signal  is  input  to 
the  low  frequency  assembly  on  multimode  fiber.  The  fiber-optic  receiver  output  is  ±  15  kHz 
bandpass  filtered  to  eliminate  aliasing  of  spectral  components,  then  routed  to  the  output  matrix 
switches.  One  output  of  the  matrix  switch  feeds  a  zero  crossing  detector  for  1-second  phase 
measurements.  The  zero  crossing  detector  generates  a  1  PPS  output  that  is  routed  over  fiber 
optics  to  the  frequency  counter.  The  other  outputs  of  the  matrix  switch  are  sent  on  coax 
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488  bus  is  used  to  communicate  with  two  frequency  synthesizers  and  a  time  interval  counter. 
(One  synthesizer  supplies  the  local  oscillator  for  the  last  downconversion  to  1  Hz  as  shown 
in  Figure  4,  while  another  supplies  the  sample  clock  for  A-D  conversion.)  The  VME  chassis 
contains  a  Skybolt  8116-V  vector  processor  and  an  Analogic  DVX  2503  16-bit,  400  kHz  A-D 
converter. 

The  Skybolt  computer  is  delivered  with  its  own  Unix-based  operating  system,  which  allows  the 
execution  of  one  user  program.  We  have  written  the  one  user  program  to  provide  custom 
real-time  multitasking  and  digital  signal  processing.  The  program  is  designed  to  accomodate 
one  test  at  a  time,  in  the  form  of  an  execution  script  including  the  digital  signal  processing, 
along  with  some  small  Skybolt  system  tasks.  The  code  is  written  in  C  and  Fortran. 

The  software  on  the  Sun  runs  with  the  Unix  operating  system  using  a  Motif-style  window 
manager  environment.  Custom  screens  allow  operators  to  use  the  stability  analyzer  with  only 
occasional  reference  to  an  instruction  manual.  Unique  test  script  files  are  compiled  at  run-time 
to  control  test  tasks,  which  are  started  in  the  Sun  and  executed  in  the  Skybolt.  The  scripts 
are  written  in  a  custom  language,  similar  to  Structured  Query  Language  (SQL),  including 
higher-level  operations  such  as  Define,  DoWhile,  If,  etc.  The  Sun  code  is  written  in  C,  some 
of  which  is  computer  generated  by  programming  tools  and  utilities,  mainly  Builder  Xcessory, 
Lex,  and  Yacc. 

The  signal  processing  software  supports  tests  for  Allan  deviation  of  phase  and  differential  phase, 
time  series  of  phase  and  amplitude,  and  spectra  of  signal,  phase,  and  amplitude.  Each  of  17 
distinct  tests  can  be  selected  by  the  operator  with  a  single  mouse  click  on  the  display.  The  test 
configuration  parameters  (input  source,  sample  rate,  averaging  time,  etc.)  are  automatically 
loaded  from  editable  configuration  files,  and  can  also  be  modified  at  the  display  by  the  operator. 

The  sample  clock  for  A-D  conversion  comes  from  a  Hewlett  Packard  3325A  synthesizer,  referred 
to  10  MHz  from  the  Reference  Frequency  Distribution  Assembly.  Although  the  A-D  converter 
can  handle  400  kHz,  the  limit  of  the  current  implementation  is  230  kHz.  Nevertheless,  this 
rate  is  adequate  to  handle  two  of  the  widest  baseband  signals  (bandwidth  45  kHz)  from  the 
Deep  Space  Network  Radio  Science  open-loop  receiver.  The  frequency  span  of  spectra  can 
vary  from  50%  of  the  sample  rate  down  to  an  arbitrarily  small  band  about  the  carrier. 

Digital  Signal  Processing  (DSP)  Algorithms 

Vectorized  Processing 

The  signal  processing  routines  run  on  a  single-board  computer,  the  40  MHz  Skybolt,  containing 
an  Intel  I860,  a  floating-point  vector  processor  with  its  own  high-speed  data  cache.  To  achieve 
the  best  computational  throughput  on  this  processor,  we  avoided  recursive  operations,  such  as 
phase-locked  loops  and  recursive  digital  filters,  in  favor  of  sequential,  nonrecursive  operations 
on  large  arrays,  such  as  element-by-element  vector  arithmetic,  inner  products,  finite-impulse- 
response  (FIR)  digital  filters,  and  the  fast  Fourier  transform  (FFT),  all  of  which  are  supported 
by  Sky  Computer’s  vector  library  and  compiler.  Throughputs  of  25-30  million  floating-point 
operations  per  second  were  achieved. 
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cable  to  A-D  converters  in  the  VME  Assembly  for  other  measurements  of  signal,  phase,  and 
amplitude.  A  digital  I/O  assembly  receives  RS232  switch  commands  from  the  computer  to 
address  the  switch  decoders  that  actuate  the  matrix  switches  and  set  attenuation  values. 

Zero  Crossing  Detector 

The  stability  analyzer  employs  two  methods  for  phase  detection:  one  method  using  software 
processing  of  A-D  samples,  and  another  using  a  time  interval  counterisi.  For  the  second 
method,  we  use  a  new  design  of  zero  crossing  detector  that  has  reduced  time  jitter  compared 
to  previous  designs^l.  In  operation,  the  zero  crossing  detector  heterodynes  the  signal  to  1  Hz, 
then  processes  the  1  Hz  output  to  produce  1  Hz  rate,  30  microsecond-wide  pulses  that  are 
sent  over  fiber  optics  to  the  time  interval  counter. 

Time  Interval  Counter 

A  HP  5334B  Counter  is  modified  to  accept  inputs  from  the  rear  panel,  and  to  accept  the 
fiber-optic  .signal  from  the  zero  crossing  detector  and  a  10  PPS  signal  from  the  reference 
distribution  assembly. 

Reference  Frequency  Distribution 

This  assembly  distributes  a  high-stability  10  MHz  station  reference  to  the  frequency  synthesizers 
and  the  time  interval  counter,  and  also  generates  a  10  pulse  per  second  signal  used  by  the  time 
interval  counter  for  phase  detection. 

Environmental  Concerns 

The  stability  analyzer  has  been  designed  to  minimize  influence  of  the  environment  on  mea¬ 
surements.  The  most  environmentally  sensitive  equipment  is  placed  in  the  frequency  standards 
room  where  ambient  temperature  stability  is  better  than  ±  O.FC.  All  signals  between  the  100 
MHz  assembly  and  the  stability  analyzer  racks  are  connected  through  fiber  optics  to  eliminate 
groundloops  that  could  induce  powerline  spurious  into  measured  output.  The  analog  electronics 
of  the  100  MHz  and  low  frequency  assemblies  are  temperature  stabilized  with  a  thermoelectric 
control  sy.stem  that  reduces  room  temperature  variations  by  a  factor  of  20.  Magnetic  shields 
around  the  electronics  attenuate  magnetic  fields  by  more  than  20  dB,  thereby  minimizing  pickup 
of  AC  powerline  harmonics. 


Controller  and  Software  Design 

The  Controller  Assembly  consists  of  a  Sun  Microsystem  Sparc  2  general  purpose  computing 
.system,  with  an  attached  VME  computer  chassis.  The  Sparc  2  performs  the  user  interface 
function,  hardware  control,  and  the  display  and  logging  of  test  results.  The  computer  includes 
an  Integrix  SBus  expansion  unit,  a  1.2  GByte  hard  disk,  a  5.0  Gbyte  Exabyte  tape  drive,  a  CD 
ROM  reader,  along  with  the  usual  monitor,  keyboard,  mouse,  and  laser  printer.  Serial  ports 
are  used  for  communication  with  the  analog  hardware  and  a  time  code  translator,  and  an  IEEE 
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Sampling  the  Video  Signal 

We  discuss  here  only  the  processing  of  the  signal  through  the  A-D  converter;  the  processing 
of  1-Hz  zero-crossing  signals  through  the  counter  has  previously  been  documentedisi .  The 
analog  “video”  signal  is  specified  to  be  a  sinewave  with  weak  sidebands  in  a  known  frequency 
band  about  the  carrier.  (The  total  sideband  power  should  not  exceed  about  —30  dBc.)  First, 
this  signal  has  to  be  sampled  at  a  such  a  rate  that  the  sidebands  of  the  digitized  signal 
faithfully  reproduce  the  sidebands  of  the  analog  signal.  For  example,  the  output  of  the  100 
MHz  Interface  Assembly  is  a  100  kHz  signal  with  sidebands  between  85  kHz  and  115  kHz. 
If  this  is  sampled  at  80  kHz,  the  sampled  signal,  which  lives  in  a  40  kHz  band,  has  a  carrier 
at  20  kHz  and  sidebands  between  5  and  35  kHz.  The  16-bit  A-D  necessarily  adds  its  own 
noi.se  and  distortion;  fortunately,  by  adjusting  the  sample  rate  one  can  reduce  their  effects  on 
measurement  results  by  whitening  the  noise  and  moving  the  aliased  harmonic  distortion  images 
away  from  the  frequency  band  of  interest. 

Overview  of  Signal  Processing 

To  allow  the  user  to  check  the  overall  quality  of  the  signal,  we  supply  a  test  called  “full 
band  spectrum”.  This  test  simply  computes  a  spectrum  of  the  sampled  signal  in  the  maximum 
frequency  span  available,  namely,  half  the  sample  rate  /^.  Also  provided  are  snapshot  plots  of 
A-D  samples  vs  time. 

The  main  job  of  the  DSP  is  to  extract  the  phase  and  amplitude  modulations  from  the  digitized 
video  signal  within  a  user-selected  frequency  B  of  the  carrier.  Two  processes  for  this  are 
supplied,  called  medium  band  and  narrow  band.  Medium  band  is  used  for  B  from  fs/4  down 
to  /s/256.  Narrow  band  is  used  for  for  smaller  values  of  B,  with  essentially  no  lower  bound 
except  that  implied  by  the  user’s  patience.  These  processes  are  de.scribed  below.  First,  however, 
we  describe  a  vectorized  algorithm  for  sinewave  analysis  that  underlies  much  of  the  processing. 

The  Pony  Computation 

At  the  heart  of  the  DSP  is  a  simple  vectorized  algorithm  for  estimating  the  frequency,  phase, 
and  amplitude  of  one  batch  of  a  sampled  sinewave.  It  was  obtained  by  adapting  Prony’s  method 
of  harmonic  analysis!*!  to  the  case  of  just  one  harmonic  component,  the  carrier  itself.  Given  an 
A’-point  data  vector  {xn,  n  =  0,. . .  V-1),  we  wish  to  fit  a  sampled  sinewave  Cn  =  A  cos(wn  +  0). 
The  computation  is  divided  into  two  parts:  Pony  1,  which  estimates  frequency  w,  and  Pony 
2,  which  estimates  A  and  0.  The  Pony  1  computation  uses  the  observation  that  the  noiseless 
sinewave  Cn  satisfies  the  difference  equation  =  (2  cos  u)cn.  Accordingly,  we  estimate 

2  cos  uj  as  the  regression  coefficient  of  the  vector  {xn-\  +2:„+i)  on  the  vector  where  n 
runs  from  1  to  A  —  2.  This  computation  requires  only  two  inner  products,  of  form  ^^id 
Y^XnXn+1,  plus  somc  Scalar  arithmetic.  For  use  in  Pony  2  and  elsewhere,  we  also  generate  a 
complex  vector  of  powers  u”,  where  u  =  exp(— iw),  by  means  of  a  vectorized  “powers”  algorithm 
that  takes  advantage  of  the  Skybolt  architecture. 

Pony  2  uses  w  to  estimate  A  and  B  by  solving  the  two-parameter  least-squares  problem 
Xn  =  a  cos  u)n  —  b  sin  un  for  the  unknowns  a  and  b.  The  only  vector  computation  needed  is 


are  combined  by  subtraction  (with  some  adjustments)  to  give  differential  phase,  which  can  be 
post-processed  in  the  same  way  as  single-channel  phase  residuals. 


Spectral  Estimation 

Direct  FFT-based  spectral  estimation  methods  are  iisedl’l.  The  sequence  of  operations  applied 
to  a  data  buffer  is  detrending,  tapering,  zero-padding  to  a  power-of-2  FFT  size,  applying  a 
real  or  complex  FFT,  squaring  the  magnitude,  equalizing  the  lowpass  decimation  filter,  and 
scaling.  Some  of  these  elements  are  discussed  below.  A  sequence  of  spectral  estimates  can  be 
averaged  to  produce  a  run  spectrum  with  greater  statistical  stability. 

Spectral  density  of  signal,  phase,  or  fractional  amplitude  deviation  is  displayed  in  units  of 
dBc/Hz,  i.e.,  single-sideband  power  per  Hz  relative  to  total  (carrier)  power,  expressed  in 
decibels.  Thus,  a  phase  spectrum  shows  £(/)  =  S^{f)/2.  A  signal  spectrum  shows  both 
sidebands. 

Each  spectrum  produced  by  the  analyzer  has  an  associated  resolution  bandwidth  b,  which  is  just 
the  two-sided  noise  bandwidth  of  the  spectral  window.  The  power  of  a  narrow  spectral  line  in 
dBc  equals  its  displayed  level  in  dBc/Hz  plus  101og5o6.  Both  6  and  lOlog^o^  are  reported  to 
the  user. 

Detrending 

Before  applying  the  FFT  to  a  data  array,  this  analyzer  preconditions  the  array  by  subtracting 
a  linear  fit  obtained  by  drawing  a  straight  line  between  the  centroids  of  the  first  sixth  and 
the  last  sixth  of  the  graph  of  data  vs  time.  This  procedure  removes  both  the  level  and  slope 
divergencies  characteristic  of  certain  processes  with  stationary  second  increments!*®!,  and  allows 
the  average  of  many  array  spectra  to  converge  to  a  stable  run  spectrum.  This  avoids  a  problem 
noticed  by  Walls,  Percival,  and  Irelanl**!,  who  preconditioned  their  data  by  subtracting  the 
mean;  they  found  that  the  estimated  spectrum  for  noise  with  a  true  f~‘^  spectrum  depended 
on  the  number  of  array  spectra  that  were  averaged.  For  full  band  spectnim,  no  detrending  is 
needed  because  most  of  the  energy  is  in  the  carrier. 

Data  Tapering 

To  avoid  problems  of  energy  leakage  from  high  portions  of  a  spectrum  into  lower  portions,  each 
data  array  is  multiplied  by  a  tapering  sequence  drawn  from  a  family  of  functions  called  discrete 
prolate  spheroidal  sequences  (DPSS).  (Actually,  we  use  a  set  of  convenient  approximations,  the 
“trig  prelates”  developed  by  GreenhalK*^!.)  For  full  band  and  medium  band  spectra,  we  use 
a  single  bell-shaped  taper  from  this  family.  For  narrow  band  spectra  we  use  a  nonadaptive, 
unweighted  version  of  Thomson’s  multiple-taper  method!*^ ’!.  An  array  of  detrended  data  is 
tapered  by  four  orthogonal  tapering  sequences,  giving  rise  to  four  distinct  “eigenspectra”,  So(/) 
through  S-slf).  These  are  averaged  to  produce  the  spectral  estimate  S(/)  for  the  array.  In  a 
broadband  noise  region,  the  Sk{f)  are  approximately  uncorrelated,  and  hence  S(/)  has  about 
one-fourth  the  variance  of  each  Sk{f).  For  a  given  frequency  resolution,  the  desired  statistical 
stability  is  achieved  from  fewer  data  arrays. 
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Then  A  and  9  are  obtained  from  Aexp(t^)  =a  +  ib. 


Medium  Band  Processing 

This  mode  of  processing  operates  by  a  sequence  of  mixing  and  filtering  to  extract  the  complex¬ 
valued  analytic  signal,  containing  only  the  power  from  the  positive-frequency  side  of  the 
original  waveform  PI,  from  which  the  amplitude  and  phase  modulations  can  be  extracted  by  a 
rectangular-to-polar  operation.  The  Pony  1  calculation  estimates  the  carrier  frequency  f^,  and 
a  mixing  signal  exp{—i2'Kfct)  is  generated  by  the  powers  algorithm.  After  the  right-hand  part  of 
the  carrier  is  mixed  to  zero  frequency,  a  FIR  lowpass  decimation  filter  is  applied  to  eliminate 
the  other  part  of  the  carrier  and  to  select  the  desired  frequency  span  (-B,  B).  The  result  is 
the  desired  analytic  signal  within  B  of  the  carrier,  shifted  to  zero  frequency. 

The  analytic  signal  is  the  basis  of  all  further  processing.  If  a  signal  spectrum  is  wanted,  then 
a  two-sided  spectnim  is  generated  after  removing  the  DC  component  (the  shifted  carrier).  If 
amplitude  or  phase  are  wanted,  then  a  rectangular-to-polar  operation  is  applied  and  the  phase 
sequence  unwrapped  from  {-it,  tt).  (Wan,  Austin,  and  Vilarf*!  give  a  more  efficient  unwrapping 
method.) 

Narrow  Band  Processing 

In  this  mode  of  computation,  the  A-D  data  for  the  whole  run  are  processed  in  contiguous 
batches  of  size  N,  each  of  which  is  analyzed  by  both  parts  of  the  Pony  computation  to  produce 
a  sample  of  batch-averaged  frequency,  amplitude,  and  phase.  The  bandwidth  of  the  extracted 
amplitude  and  phase  samples  is  fs/{2N).  Because  of  the  efficiency  of  the  Pony  computation, 
the  DSP  can  keep  up  with  the  stream  of  A-D  samples  at  the  highest  rate  of  the  A-D  converter, 
400  kHz,  although,  as  mentioned  above,  the  analyzer  is  currently  limited  to  a  total  sample  rate 
of  230  kHz. 

For  computational  efficiency',  N  has  to  be  at  least  200.  To  save  storage,  we  allow  batches  no 
greater  than  a  designated  maximum  batch  size  (now  8192).  Because  we  also  wish  to  allow 
arbitrarily  small  analysis  bandwidths,  the  batch  averages  can  themselves  be  averaged  together  in 
groups  of  arbitrary  size  r  to  produce  samples  with  bandwidth  f,,/{2Nr).  In  choosing  this  crude 
lowpass  decimation  method,  we  accepted  some  aliasing  problems  to  gain  simplicity,  consistency, 
and  efficiency. 

The  phase  information  computed  by  Pony  1  and  2  is  local  to  each  batch,  and  is  known  modulo 
2ir  only.  We  have  devised  an  algorithm  to  process  these  local  data  into  a  sequence  of  global 
phase  residuals;  it  is  essentially  the  same  as  the  algorithm  used  for  processing  the  1  Hz  zero 
crossing  counter  readingsl®!.  For  the  algorithm  to  succeed,  the  frequency  must  be  changing 
slowly  enough  from  batch  to  batch  so  that  the  current  batch  phase  can  be  predicted  from 
earlier  ones  within  tt.  The  algorithm  issues  an  alarm  if  any  prediction  error  exceeds  7r/2  in 
absolute  value. 

The  low-rate  sequence  of  amplitude  and  phase  residuals  extracted  by  the  narrow  band  process 
can  be  subjected  to  a  variety  of  post-processing  functions,  including  time-series  display,  spectral 
estimation,  and  Allan  deviation.  For  a  two-channel  test,  the  phase  residuals  of  the  two  channels 
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Allan  Deviation 


From  a  stream  of  narrow  band  or  1  Hz  zero  crossing  phase  residuals,  the  analyzer  produces 
estimates  of  Allan  deviation  with  estimated  drift  removed,  using  the  simple  three-point  drift 
estimator  recommended  by  Weiss  and  Hackmanl^'^l.  The  required  r-overlapped  sums  for  first 
and  second  moments  of  second  r-differences  of  phase  are  accumulated  in  real  time. 

To  generate  conservative  error  bars  for  plus  or  minus  one  standard  deviation  of  Allan  variance, 
we  assumed  a  random-waik-freqiiency  model  of  phase  noise.  Using  a  method  of  Greenhallf^^l, 
we  carried  out  a  numerical  computation  of  u,  the  equivalent  degrees  of  freedom  of  the  drift- 
removed  Allan  variance  estimator,  as  a  function  of  M,  the  number  of  summands.  The  sequence 
of  (/  vs  M  was  fit  with  a  simple  empirical  formula.  Then,  if  a  is  the  estimated  Allan  deviation, 
the  reported  error  bar  is 


Because  of  severe  negative  bias  of  the  drift-removed  estimator  for  .small  M,  results  are  reported 
only  for  M  >  4. 


Test  Methods 

A  series  of  tests  of  the  stability  analyzer  were  conducted  at  JPEs  Frequency  Standards  Laboratory 
in  order  to  demonstrate  first,  that  the  results  of  the  stability  analyzer  agree  with  those  of  other 
measurement  equipment,  and  second,  that  it  meets  its  noise  floor  requirements.  Noise  floor 
results  are  given  in  Table  1. 

Allan  deviation  runs  of  at  least  24  hours  duration  were  carried  out  on  pairs  of  100  MHz 
frequency  standards.  The  results  were  compared  to  those  from  an  existing  FSL  Allan  Deviation 
test  set  and  found  to  agree  within  5%.  The  noise  floor  was  measured  by  splitting  the  single 
output  of  an  H-maser  and  applying  it  to  two  inputs  of  the  stability  analyzer.  These  tests  were 
carried  out  in  both  zero  crossing  detector  mode  and  the  narrow  band  phase  modes. 

Time  series  of  differential  phase  were  tested  using  a  HP  3326  dual  channel  synthesizer  as  the 
input  source.  The  two  outputs  of  the  synthesizer  were  manually  steered  in  frequency  to  produce 
phase  drifts  of  known  amplitude.  Comparison  was  made  to  the  results  from  a  HP  8508  phase 
meter,  and  also  to  a  strip  chart  recording  the  phase  difference.  This  last  signal  was  developed 
by  simple  mixing  between  the  two  outputs  of  the  synthesizer.  These  tests  were  also  run  with 
both  channels  of  the  synthesizer  set  at  the  same  frequency  for  at  least  15  hours,  to  observe  the 
noise  floor. 
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Table  1.  Stability  Analyzer  Noise  Floors 


Test 

Input  Source 

Allan  Deviation 

100  MHz 

tau 

sigma 

1  sec 

10  sec 

100  sec 

1000  sec 

6  X  lO'^'" 

2  X  10-1''^ 

2  X  10-1“ 

3  X  10- ‘i" 

Phase  Spectrum 

100  MHz 

Freq. 

Spectral  Density 

1  Hz 

10  Hz 
>100  Hz 

-126  dBc^Hz 
-135  dBc/Hz 
-142  dBc/Hz 

Signal  Spectrum 

baseband 

1  Hz 

-92  dBc/Hz 

10  Hz 
>100  Hz 

-97  dBc/Hz 
-98  dBc/Hz 

Phase  Spectrum 

1  Hz 

-98  dBc/Hz 

10  Hz 
>100  Hz 

-104  dBc/Hz 
-105  dBc/Hz 

Amplitude  Spectrum 

1  Hz 

-70  dBc/Hz 

10  Hz 
>100  Hz 

-85  dBc/Hz 
-88  dBc/Hz 

Diff.  Phase 

baseband 

Avg,  Time 

Phase  Error 

1  sec 

1000  sec 

<0.001  deg  rms 
<0.04  deg  rms 

Spectra  were  tested  in  a  variety  of  ways.  The  signal  sources  were  two  H-masers  for  the  100 
MHz  inputs,  one  H-maser  and  an  HP  8662  synthesizer,  or  one  or  two  HP  3325  synthesizers 
for  the  ba.seband  analog  inputs,  In  the  latter  two  cases,  one  synthesizer  was  modulated  either 
by  another  synthesizer  to  simulate  spurious  signals,  or  by  a  HP  3561  noise  source  to  simulate 
phase  noise.  The  spectrum  was  then  compared  to  the  results  from  a  HP  3589  or  3561  spectrum 
analyzer.  The  results  agreed  within  a  typical  2  dB  peak-to-peak  variation  between  spectral 
bins.  For  noise  floor  tests,  a  single  H-maser  signal  was  divided  and  applied  for  comparison  at 
two  inputs. 
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Fig.  1.  Stability  analyzer  in  a  typical  DSN  installation. 


Fig.  2.  Stability  analyzer  rack  arrangement. 
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